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We study the Landau damping in Bose-Fermi superfluid mixture at finite temperature. We find 
that at low temperature, the Landau damping rate will be exponentially suppressed at both the 
BCS side and the BEC side of Fermi superfluid. The momentum dependence of the damping rate 
is obtained, and it is quite different from the BCS side to the BEC side. The relations between our 
result and collective mode experiment in the recently realized Bose-Fermi superfluid mixture are 
also discussed. 


I. INTRODUCTION 

The quasiparticle is an important concept in mod¬ 
ern many-body physics. The statistics and dispersion 
of quasiparticles determine various low energy properties 
of the quantum many-body system. In the general situa¬ 
tion, there exsits interaction between the quasiparticles. 
That leads two major effects. One is the modification of 
the dispersion of excitations, corresponding to the real 
part of the self-energy of quasiparticle. Second, the in¬ 
teraction can also damp a given excitation, i.e. giving 
quasiparticle a finite lifetime, which can be reflected by 
the imaginary part of the self-energy. The damping of 
low energy excitations is responsible for many interest¬ 
ing phenomena of many-body system, such as transport 
and tliermalization. 

For example, in a uniform condensed Bose gas with 
short-range interaction, the low energy excitation is 
phonon-like quasiparticle with linear dispersion when 
wavelength is larger than the healing length. While when 
wavelength is smaller than healing length, the excitation 
is free-particle-like quasiparticle with quadratic disper¬ 
sion. At zero temperature, due to the residual interac¬ 
tion between excitations, a quasiparticle in the BEC can 
be damped by decaying into two quasiparticles. Since 
each product of decay must have an energy lower than 
the original quasiparticle, the final state phase space is 
restricted. This gives the damping rate a very sensitive 
dependence on the initial momentum 7 ~ A ; 5 [Ti. This 
decay process is the so-called Beliaev damping. At finite 
temperature, a given quasiparticle can also be damped by 
absorbing thermal quasiparticles. As a result, it is very 
sensitive to temperature, 7 ~ T 4 fc (?j. This mechanism is 
known as Landau damping, which was first discussed in 
plasma oscillation by Landau [31. It plays an important 
role in various phenomena such as anomalous skin effect 
in metals and the damping of phonons in solids. Landau- 
Beliaev damping in a uniform Bose superfluid has been 
widely studied both theoretically 00112 ] and exper¬ 
imentally |S]. In trapped system, since the low energy 
excitations are discrete, Beliaev damping is forbidden. 
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The damping of low energy modes is attributed to Lan¬ 
dau process, and the experiment damping rate has been 
found to be consistent with the theory of Landau damp¬ 
ing mm- Landau-Beliaev damping has also been stud¬ 
ied in dipolar BECs El and in the mixture of BEC with 
normal Fermi gas [12] . 

Recently a Bose-Fermi superfluid mixture has been 
first realized by ENS group m The dipole mode of this 
new superfluid mixture has been measured. It exhibits a 
frequency shift and an unusual damping behavior. This 
experimental development triggers many investigations 
on Bose-Fermi superfluid mixture nu ns] [m nz] iiH] mi 

The quasiparticles in this Bose-Fermi superfluid mix¬ 
ture have a quite unique feature. There are two gapless 
bosonic modes, which are Goldstone modes of Bose and 
Fermi superfluids. It also has a gapped fermionic mode, 
corresponding to the Cooper pair breaking. Moreover, in 
the ENS experiment, the Fermi superfluid can be tuned 
from the BCS side to the BEC side by Feshbach reso¬ 
nance. During the crossover, the behavior of three kinds 
of excitations gradually changes from the BCS limit to 
the BEC limit. In the BCS limit, the velocity of Gold- 
stone mode in Fermi superfluid is quite large, approach¬ 
ing VF/y/3 ED], while the gap of the fermionic mode is ex¬ 
ponentially small. When it is tuned to the BEC side, the 
velocity of the Goldstone mode decreases monotonously 
mm, while the gap becomes larger and larger. The 
dispersions of the three excitations at both sides are plot¬ 
ted in Fig [I] 

To understand the unusual damping behavior of dipole 
mode in the ENS experiment, Zheng and Zhai study the 
Beliaev damping of bosonic mode in Bose superfluid by 
considering its interacting with quasiparticles in Fermi 
superfluid at zero temperature El- They found that 
Beliaev process will be activated only if the excitation 
momentum exceeding a critical value k c . This threshold 
damping behavior, i.e. 7 ~ (k — k c ) a , is quite different 
at the BCS side and the BEC side. To be specific, at the 
BCS side a = 0, while at the BEC side a = 3. This is 
because at the BCS side, the damping is dominated by 
decaying into fermionic pair-breaking modes in Fermi su¬ 
perfluid. The final state phase space of fermionic mode is 
restricted by density-of-state near the Fermi surface, so 
that the damping rate is nearly a constant. On the other 
hand, such a restriction does not exist at the BEC side, 
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where the damping is dominated by decaying into Gold- 
stone modes in Fermi superfluid. So the damping rate 
grows rapidly with momentum. However, they only con¬ 
sider the zero-temperature case, in which Landau damp¬ 
ing is frozen. It is nature to ask the question, what will 
happen at finite temperature when Landau damping is 
activated. 

In this paper, we investigate Landau damping of 
bosonic quasiparticles in Bose superfluid due to its in¬ 
teracting with quasiparticles in Fermi superfluid at finite 
temperature. We consider the typical cold atom situ¬ 
ation, in which Fermi superfluid is in the strongly in¬ 
teracting regime, while Bose superfluid is in the weakly 
interacting regime. We find that unlike Beliaev damping, 
the critical momentum for Landau damping is zero, i.e. 
k c = 0 , both at the BCS side and the BEC side. We 
obtained the temperature dependence of Landau damp¬ 
ing, showing that the damping rate is exponentially sup¬ 
pressed at temperature low enough at both sides. This 
is indeed a departure from the T 4 dependence in the 
single component BEC. The momentum dependence of 
the damping rate is also obtained, and its behavior is 
quite different from the BCS side to the BEC side. At 
BCS side, the damping rate grows linearly with momen¬ 
tum, while at the BEC side, it is nearly a constant. This 
damping behavior revealing different dominated low en¬ 
ergy quasiparticles at each side. 

This paper is organized as follows: In Sec. II we con¬ 
struct the model for Bose-Fermi superfluid mixture and 
its mean field treatment. In Sec. Ill Landau damping 
is considered by perturbation method at the BCS side. 
In Sec. IV Landau damping is calculated at the BEC 
side. In Sec. V, our main results are summarized and 
the connection to experiment is discussed. 


II. THE MODEL 


a magnetic field is near a Feshbach resonance between 
fermions, gb and gbf are generally in the weakly inter¬ 
acting regime and are approximately constant according 
to the experiment setup. Therefore for Hb , we take the 
Bogoliubov mean-field theory. For Hf, we take the BCS- 
BEC crossover mean-field theory. After mean-field treat¬ 
ment, one obtains 

ij b MF = E^4« k, ( 2 ) 

k 

Hf F = J2 E k(/pk + 7l7k), (3) 

k 


where E £ = y (e£. + 2 g b n b ), d*, is the quasiparticle 
operator for bosonic Goldstone mode in Bose superfluid. 

E ^ = yj (e£ — Hf ) 2 + A 2 . / 3 k and y k are the quasiparti¬ 
cle operators for fermionic pair-breaking mode in Fermi 
superfluid. e k = h 2 k 2 /(2m*), i = b, f denotes the ki¬ 
netic energy of bosons and fermions. The corresponding 
Bogoliubov transformations of these operators are given 
by 


6k 

b * b 

= !i k a k - V k al k , 

(4) 

Ck,f 

= + 

(5) 

Ck,4. 

= u L^~ v L^-k- 

(6) 


Here momentum-dependent coefficients w k and w k are 
given by 


<K) = \h 






■4 - aA 
eL \ 


(7) 

( 8 ) 


Consider a homogeneous mixture of Bose and Fermi 
superfluid. The Hamiltonian of the superfluid mixture 
has three parts: H = Hb + Hf + Hbf, 

H b = J d 3 r |6 t (r)iL 0jb S(r) + y fe t (r)& t (r)fe(r)S(r)|, 

H f = I d 3 r |^ 4 (r)-ffo,/c CT (r) 

( r )c| ( r )c4,(r)c T (r)} , 

Hbf = 9 bf d?rtf(r)b(r)ct(r)c a (r), (1) 

a ** 

where H 0 ,i = — h 2 '\ 7 2 /( 2 m i ) — Hi, i = b, f denotes 
bosons and fermions, a =t,4- denotes the spin com¬ 
ponents of fermions. The coupling constants are re¬ 
lated to the scattering lengths: 1 /gb = mb / (Anh 2 ab) and 
1/5/ = fnf / (Anh 2 a /) + ^ k l/(2e£), where the scatter¬ 
ing lengths can be tuned by Feshbach resonance. When 


For fermionic pair-breaking mode in Fermi superfluid, 
when —1/kpCLf gets smaller from the BCS side to the 
BEC side, A will increase and \Xf will decrease. Apart 
from the pair-breaking mode, there is also a bosonic mode 
of the center-of-mass motion of Cooper pairs in Fermi 
superfluid, which is beyond BCS-BEC mean-field theory. 
Its dispersion relation at low energy is linear: Eg 1 = 
hcjk. From the BCS side to the BEC side, c/ evolves 

from Vf/V 3 to Tth 2 a m n m /m 2 [2T] |22| . where a m = 

0.6a/ | 23 J, and n m = n/ = nj. = nf. 

Therefore there are three different excitations in Bose- 
Fermi superfluid mixture: the bosonic Goldstone mode 
in Bose superfluid, whose dispersion is given by E £, the 
fermionic pair-breaking mode in Fermi superfluid, whose 
dispersion is given by E ^ and the bosonic Goldstone 
mode in Fermi superfluid, whose dispersion is given by 
Ey?. At the BCS side, the fermions form Cooper pairs 
and become superfluid of BCS type. At the BEC side, 
the fermions form strongly bound molecules and become 
superfluid of BEC type. The dispersions of the three 
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FIG. 1: Schematic of dispersions of the bosonic Goldstone 
mode in Bose superfluid (E ^), the fermionic pair-breaking 
mode in Fermi superfluid(i5 k ) and the bosonic Goldstone 
mode in Fermi superfluid ( E (a) is in the BCS side and 
(b) is in the BEC side. 


excitations at different sides are shown respectively in 
Fig. 0a) and (b). In the experiment setup, the Bose 
gas is so dilute that we can take the free-particle limit, 
i.e. E £ « £ k for the Bogoliubov mode in the Bose su¬ 
perfluid. In this case, Landau-Beliaev damping of the 
bosonic mode in Bose superfluid due to its interaction 
with itself can be ignored. 


III. DAMPING AT THE BCS SIDE 


Consider the interaction between bosons and fermions. 
Due to the existence of boson condensate, bo and can 


be treated as c-numbers, i.e. bo = b^ = \/Nb. 


the Hamiltonian by the order of '/Nt ,, we have 
H$+H$ + H§\ where 


Expand 
H b f = 


Hbf ~ 9bf Tib ^k,cT^k,c7, 


f>( 2 ) _ dbfVM 

M bf - 


V 


X] (Cq+k.o-Cq.cr^k + ll-C.), 


k^ 0 ,q ,a 


o-(3) _ ffbf_ V- -f A _ 
n bf ~ y C p-q,a C p-k,CT0 q C>k- 


k,q7^0,p,cr 


(9) 

(10) 

( 11 ) 


The leading term Hjf will only shift the chemical poten- 
tial of fermions. The subleading term will induce 

the damping of boson quasiparticles. The last term 
is a two particle scattering process, which is less irnpor- 
tant compared to H^j , so that it can be ignored. We get 

several damping channels by expressing with mean- 

~ (2) A 

field quasiparticle operators: H^f ss H\ + H- 2 + H 3 . Here 


Hi 

Hi 

H 3 


9bfVNb 

V 


x« 


V k)( U q+k U q - W q+kG^q+k^q«k + h.C. 


gbfVNb 

V 




9bfy/Nb 

V 


E(4- 


V k)(«q+k U q - W q+k i; q)7q+k7q«k + h.C. 

q^q - W k-q U q)7q/3Lq d 'k + h.C. 


( 12 ) 

(13) 

(14) 


We have ignored terms like /3 q /3 q -kdk that violate the 
conservation of energy, which will not contribute to 
damping process. 

We note that H 3 is the decay process of bosonic mode 
in Bose superfluid, which is Beliaev damping. H i and 
H 2 are the scattering of the bosonic mode in Bose su¬ 
perfluid by thermal excitations, corresponding to Landau 
damping. Unlike Landau-Beliaev damping in single com¬ 
ponent Bose gas, there exists nonzero damping threshold 
for some of the damping channels. While due to the 
energy gap A in /3 q +k and 7 q +k- The damping thresh¬ 
old for H 3 is 0 3 (k) = min[£T + U k _ q ] > 2A which is 
nonzero. That is to say, for low energy excitations with 
E £ < fl 3 (k), there will be no damping contributed from 


channel H 3 . So the Beliaev damping has a critical mo¬ 
mentum. The damping threshold of Hi and H 2 are the 
same, which is fii(k) = fl 2 (k) = min[.E q+k — Ef\ = 0. So 
the Landau damping we considered here has no damping 
threshold, i.e. the critical momentum for the damping is 
zero. 

As shown in Fig. 0a), there is also a phonon-like 
bosonic mode of center-of-mass motion of Cooper pairs 
in Fermi superfluid. We will discuss its contribution to 
the Landau damping later. 

Since the interaction between bosons and fermions is 
weak compared to the excitation energy in both bosons 
and fermions, so that we can treat it perturbatively. Ac¬ 
cording to Fermi’s Golden Rule, the rate for the Landau 
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FIG. 2: Damping rates 7 in units of 70 (70 = Ef/K) as a 
function of k/kp with 1 /(kpCLf) = —0.5 in the BCS side, 
(a) Damping rates as a function of momentum; (b) Damping 
rates as a function of temperature. 


damping described by H\ and H 2 is given by 


7fl0 = yE l^qkl [f{E{) - /(^ +k ) 

q 

(15) 

where f(E^) = [exp (E^/T) + 1] 1 is the Fermi-Dirac 
distribution function. The matrix element is 


Afqk = 


2g bf VNb l ^ b 


V 


«-<)(' 


f f 
,uL - v J „ 


q+k'A u q+k L 'q 


). (16) 


The damping rate given by © is numerically calculated 
and the result is plotted in Fig. [2j 

In Fig. [^a), the damping rate is zero for k = 0. It is 
linear in k when k is small. The linear behavior is the 
result of matrix element and density-of-state of quasipar¬ 
ticles together. The slope is larger at higher temperature 
due to more thermal excitations. 

What is more interesting is the temperature depen¬ 
dence of Landau damping rate. As shown in Fig. j2jb) , 
when T/Tp -C 1, the damping rate shows a e~ Tp '^be- 
havior. This exponential decay is very different from the 
power law behavior of Landau damping in dilute Bose 
gas. Similar as the damping threshold discussed before, 
the conservation of energy must be satisfied in the damp¬ 
ing process. Due to the existence of the gap A in (3^ 


and 7 k, in order for quasiparticle to be damped, it must 
absorb a thermal excitation with E ^ > A. However, 
as temperature tends to zero, the distribution function 
f{Ef) tends to be e _A ^ fcsT - ) . The number of thermal 
excitations is exponentially suppressed at low tempera¬ 
ture kpT < A. Therefore the damping rate is also expo¬ 
nentially suppressed. The size of the suppressed region is 
proportional to gap. The region will be larger if the BCS 
gap is tuned to be larger. For usual Landau damping in 
dilute Bose gas, thermal excitations are gapless. So even 
at low temperature, a number of thermal excitations can 
be excited and contribute to damping, leading to a T 4 
power law temperature dependence. 


IV. DAMPING AT THE BEC SIDE 

According to the result at the BCS side, the damping 
rate contributed by the pair-breaking channel will have 
a exponentially suppressed region. In the BEC side the 
molecule is tightly bounded and the pair-breaking energy 
A is quite high. So this suppressed region is large. This 
damping channel can be neglected. 

Now the damping channel of center-of-mass motion of 
Cooper pairs in Fermi superfluid is important. A com¬ 
prehensive description of Goldstone mode in Fermi super¬ 
fluid and it coupling to Bose superfluid can be obtained 
from fluctuation theory of Fermi superfluid |24j . Here to 
expose physics in a simple way, we treat Fermi superfluid 
at the BEC side as a molecular BEC. So the Hamiltonian 
is given by three parts: H = H b + H m + H bm . H b is the 
same as before. 

H m = J d 3 r {r)H 0 m d(r) + ^d t (r)d t (r)d(r)J(r)|, 

'Hbm = 9 bm J d 3 rP(r)d(r)b\r)b(r), (17) 

where E 0 , m = -fi 2 V 2 /(2?n m ) - n m . 1 /g m = 

m m / (47r/i 2 a m ). m m = 2 irif and a m = 0.6 a/. For H m , 
we take the Bogoliubov mean-field theory. The result is 

^r=E^xixk, (is) 

k 

where E^ = + 2g m n m ), \k is the quasipar¬ 

ticle operator for bosonic Goldstone mode in molecu¬ 
lar BEC. e^ = H 2 k 2 / (2m m ) is the kinetic energy of 
molecules. The corresponding Bogoliubov transforma¬ 
tion of molecule operators is given by dk = u kXk — t’ k X_k- 
where momentum-dependent coefficients and c k l are 
given by 


<«) = ( £ * + E 9 S nm ±i)- ( 19 ) 

Due to the existence of condensates, 60 , &o, do and dj 
can be replaced by b 0 = b^ = y/Nb, do = dj = y/N m . To 
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the order of i/N b and \ r N mi H bm = H^ + + H ^, 

where 

&bm =9bm Z(nmblb k + Hfcd^dk), (20) 

k^O 

M2) _gbmVN m (A i it At t. 

Hbm — y Z_^ (^q + ^-q)^k+q^k 

k,q#0 


. Ohm y/Nb 
V 

E ( d q+k d q^k + ^q_k d q^k)’ 

(21) 



k.q^O 


Ohm 

~v~ 

E 

d q-k d q^p+k^P- 

(22) 


p.q^O.k^q 


The leading term H b m only modify the chemical po- 
tential. The subleading term H tmi will contribute to the 

damping. The last term H b j is a two particle scatter¬ 
ing process, which is less important compared to the 
and can be ignored. Express H^ with mean- 
field quasiparticle operators, we get damping channels 
H ( h f r { « Hi + Hi + H 3 + Hi, where 


Hi 


Hi 


H :i 


H 4 


Ohm y/Nm Y'V m mi, b b , b 6W , 1 

- v - q “ u q )K-q W k + <-q<)Xq<- q a k + h.C., 


V 

9brn yM 

V 


k.q 


J^( U k - t'kK-qt'qXU-q^ + h - C ‘’ 


k.q 


Obm^/' Nj 


■E(‘ 


m \^/ m _ r 


y / ,v^q+k "q+k 
k,q 


)Wk«qXq +k Q ; q d 'k + h -C-, 


E(^k - 4)« + k< + < + k<)xUxq d k + h.c„ 

k.q 


(23) 

(24) 

(25) 

(26) 


We have ignored terms like x_ q d q +kQ'k that violates the 
conservation of energy that will not contribute to damp¬ 
ing process. 

Here Hi is the decay of bosonic mode in Bose super¬ 
fluid, which is Beliaev damping. H 3 and H 4 are scat¬ 
tering by thermal excitations, which are Landau damp¬ 
ing. Special attention should be paid on H\: Xq d k_ q d k 
gives Beliaev damping, while its Hermitian conjugation 
d q+kXq d k can contribute to Landau damping. Similar as 
the BCS side, dk and y k are different quasiparticles with 
different dispersions. To satisfy conservation of energy 
and momentum, the Beliaev damping channels H 1 and 
Hi have damping thresholds for low-energy excitations: 


Oi(k) 


k 2 / ( 2 m b ) 
c m k - m. b c 2 m /2 


k < in b c m , 
k > m b c m , 


0 2 (k) = c m k, 


(27) 


respectively. The Landau damping channels have no 
damping threshold, i.e. have a zero critical momentum. 


The matrix elements of Landau damping are given by: 


M Obm VNm, m mi, b b , b bl 

l.qk — y (tiq - fq )(«k-q U k + V k-q l ’k)> 

gbmVH m . m m . , u h , .. 6.6 n 


Af 3> qk = — 
My q k = 


/ m m \ / b b i b b \ 

y ( U q+k — ^q+kX^k^q + V k U q)l 

gbmy/Hb t b b . 


V 


-K-<)« + k<+<Vk^)- (28) 


As mentioned before, we only consider the free-particle- 
limit for Goldstone mode in Bose superfluid. So n k and 
v\. can be simplified as u k « 1, c k « 0. We immediately 
get that Afi.qk = M 3iq k = 0. Physically, this is because 
in the free-boson limit d-k ~ b k- Thus Hi and H 3 violate 
the conservation of number of bosons. Therefore we only 
consider damping channel U 4 . 

According to Fermi’s Golden Rule, this Landau damp¬ 
ing rate is given by 


7(k) = Y £ [/(£™_ k ) - f(E™)] , 

(29) 

where f(E™) = [exp ( E™/T ) — 1] 1 is the Bose-Einstein 
distribution function. The conservation of energy E k + 
E™_ k — E™ = 0 will give a momentum lower bound g c (k) 



6 



FIG. 3: Damping rates 7 in units of 70 (70 = Ef/K) as a 
function of k/k F with 1 /(fra/) = 0.5 for the BEC side, (a) 
Damping rates as a function of momentum; (b) Damping rates 
as a function of temperature. 


for the thermal excitations, which is 

k f _ Hk 


Qc (k) = 


1 - 


2?77/£)Cfj 


(30) 


Only thermal excitations with q > g c (k) will contribute 
to the damping. The physical reason for this lower bound 
can be understood similar as the damping threshold as 
following: for infinitesimal q, the conservation of en¬ 
ergy and momentum cannot be satisfied simultaneously, 
because Xk are quasiparticles with different dispersions 
from dk- To be more specific, Xk has a large velocity and 
has higher energy compared to dk- 

The damping rate given by (29) is numerically calcu¬ 
lated and the result is plotted in Fig. [3j In Fig. [3ja), 
the damping rate is nonzero at k = 0. At leading order, 
the damping rate decreases with momentum. This can 
be understood as the competition of matrix element and 
the density-of-state of the quasiparticles. The leading 
and the subleading term of q(k) can be calculated an¬ 
alytically by taking the asymptotic expansion of matrix 
element. The result is 


7(k) = 


9bm n b™L 

Airffinib 


2 m 


mCm 

hTp 


T 


1 + ^T'U 

2 m b T F ) 


+0(k 2 ) 

(31) 


The intercept at k = 0 is proportional to T/Tp. This is 
reasonable because when T = 0, there are no thermal ex¬ 
citations. The subleading term is always negative, which 
agrees with the numerical calculations, see Fig. [3](a) . 


The temperature dependence of the damping rate is 
plottd in Fig. §b). We note that 7 (T/Tp) is linear in 
T/Tp at high temperature. At low temperature, when 
T/Tp < 1, 7 (T/Tp) oc e~ T¥ l T . So at the BEC side, the 
temperature dependence also shows an exponentially de¬ 
cay behavior, which is similar to the BCS side. However, 
the physical origin of this exponentially suppression is 
not the gap, but the momentum lower bound. Since only 
excitations with q > q c can contribute to the damping, 
there is an energy lower bound E c = Hc m q c for ther¬ 
mal excitations to involve into Landau damping process. 
For low temperature kpT < E c , the involved thermal 
excitations are exponentially suppressed, leading to an 
exponentially suppressed damping rate. The size of this 
suppressed region is proportional to energy lower bound 
E c . According to (30), larger the k, larger the E c , hence 
larger the suppressed region. This is agreed with Fig. 
§ b ) and is also a distinct feature from the BCS side. 
For the Landau damping in single component BEC, the 
quasiparticles are also gapless. But there is only one 
kind of quasiparticle and thus no energy lower bound. 
The conservation of energy and momentum can always 
be satisfied there. As a result, there is no exponential 
suppressed region. 


Now we would like to point out that the Landau damp¬ 
ing channels via interacting with the Goldstone mode in 
Fermi superfluid also exist at the BCS side. Although the 
detail of these channels cannot be covered at mean-field 
level, they must obey the conservation of energy and mo¬ 
mentum. The velocity of Goldstone mode in Fermi super¬ 
fluid at the BCS side is about vf/\/2> and is quite large. 
So based on our previous analysis, the corresponding en¬ 
ergy lower bound hc m q c is also very large at the BCS 
side, such that the damping process via those channels 
are highly suppressed in a large temperature range. It 
can be ignored comparing to the damping channels via 
interacting with the fermionic pair breaking modes at the 
BCS side. 


At the end of this section, we should emphasize that 
our analysis of temperature dependence in Landau damp¬ 
ing is only viable at a relative low temperature, since both 
the gap of the Fermi superfluid and the condensate frac¬ 
tion in Bose superfluid will decrease with the increasing 
of temperature. A comprehensive self-consistent analysis 
of temperature dependence should take the changing of 
the gap and the condensate fraction into account. How¬ 
ever, at low temperature both the gap and the conden¬ 
sate fraction are not sensitive to the temperature, so our 
result is still reasonable in this region. 
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V. SUMMARY 

In summary, we have studied Landau damping in Bose- 
Ferrni superfluid mixture at finite temperature. Unlike 
Beliaev damping in Bose-Fermi superfluid mixture at 
zero temperature, Landau damping has no critical mo¬ 
mentum, since any quasiparticle with infinitesimal mo¬ 
mentum can be damped by absorbing a thermal quasi¬ 
particle. However, due to energy-momentum conserva¬ 
tion, thermal excitations in Fermi superfluid will be in¬ 
volved into Landau process only if their energy is larger 
than a lower bound E c . When the temperature is lower 
than E c , Landau damping rate will be exponentially sup¬ 
pressed. For fermionic quasiparticles in Fermi superfluid, 
the lower bound is determined by pair-breaking gap A. 
For bosonic excitations, the lower bound is determined 
by Hc m q c . This exponential suppression is quite different 
from the power law temperature dependence of Landau 
damping in single component BEC. The reason behind 
is that in Bose-Fermi superfluid mixture, quasiparticle in 
Bose superfluid is damped by coupling to thermal exci¬ 
tations in Fermi superfluids with totally different disper¬ 
sion. In single component BEC, quasiparticle is damped 
by interacting with itself, so that energy-momentum con¬ 
servation can always be satisfied. Therefore there will be 
no energy lower bound and no exponential suppression. 
In dipolar Bose gas, there is also an exponential sup¬ 
pression region for Landau damping m- That is due 
to its unusual quasiparticle dispersion with maxon-roton 
structure. 


In principle, both damping channels contributed from 
bosonic quasiparticles and fermionic quasiparticles in 
Fermi superfluid exist at both sides. However, at the 
BCS side, we have A <C Iic m g c , so that Landau damping 
is dominated by fermionic excitations in Fermi superfluid. 
At the BEC side, the situation is opposite, A >> hc m q c . 
So bosonic quasiparticles in Fermi superfluid dominate 
Landau process. This result is similar to the case of zero- 
temperature Beliaev damping. The different dominated 
low energy quasiparticles also lead to totally distinct mo¬ 
mentum dependence of damping rate at the BCS side and 
the BEC side. 

In the ENS experiment m, damping of dipole modes 
has a critical momentum even the experiment is done 
at finite temperature. Based on our previous analysis, 
this phenomenon can be understood as following: the 
temperature of the experiment is sufficient low that Lan¬ 
dau damping with zero critical momentum is highly sup¬ 
pressed. Then Beliaev damping with nonzero critical mo¬ 
mentum will show up m- Our prediction of momen¬ 
tum dependence of Landau damping could be observed 
in the same experiment setup at relative high tempera¬ 
ture, when Landau damping will dominate over Beliaev 
damping. No critical momentum will be observed there. 
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